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Abstract 

An exact series expansion is now known for the Principal Field Emission Elliptic Function v, in 
terms of a complementary elliptic variable l' equal to y 2 , where y is the Nordheim parameter. 
This expansion was originally found by using the algebraic manipulation package MAPLE. 
This paper presents a formal mathematical derivation. Recently, it has been discovered that 
v(l') is a particular solution of the ordinary differential equation (ODE) Z'(l — l')d 2 v/dl' 2 = nv, 
when the index n — 3/16. This ODE appears to be new in mathematical physics and elliptic- 
function theory. The paper first uses an 1876 result from Cayley to establish the boundary 
condition that dv/dl' satisfies as I' tends to zero. It then uses the method of Frobenius to 
obtain two linearly independent series solutions for the ODE, and hence derives the series 
expansion for v(l'). It is shown that terms in ha. I' are required in a mathematically correct 
solution, but fractional powers of I' are not. The form of the ODE also implies that simple 
Taylor expansion methods cannot generate good approximation formulae valid over the whole 
range < I' < 1; this conclusion may also apply to barriers of other shapes. It is hoped that 
this derivation might serve as a paradigm for the treatment of other tunnelling barrier models 
for cold field electron emission, if in any particular case an ODE can be found for which the 
tunnelling-exponent correction function is a particular solution. 



1 Introduction 

Fowler-Nordheim (FN) tunnelling [1] is electric-field-induced electron tunnelling through a roughly 
triangular barrier. When the emission barrier is strong and the penetration coefficient small, 
there is a low-temperature emission regime (including room temperature) known as "cold field 
electron emission (CFE)" . Tunnelling and CFE are processes of significant technological interest, 
in particular for the prevention of vacuum breakdown, the development of cold-cathode electron 
sources, and internal electron transfer processes in some types of electronic device. Tunnelling and 
CFE are at the heart of several multi-billion-dollar industries. 

CFE from a metal-like conduction band is described by a family of approximate equations known as 
"Fowler-Nordheim-type" (FN-type) equations. FN-type equations for the emission-current density 
J, as a function of the local work-function <j> and the local barrier field F, have the form 

J = Xa^F 2 exp(-nbc/) 3 / 2 /F) (1) 

where a and b are the First and Second FN Constants as usually defined (e.g. [2]), and A and 
are generalised correction factors. The forms of A and \i depend on the assumptions and 
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approximations made in the derivation, in particular on the shape assumed for the tunnelling 
barrier. 

For a planar emitter, a "Schottky-Nordheim" (image-rounded) barrier model [3, 4] is normally 
used. In this case, fi is given by the mathematical parameter v-p, which is a particular value (see 
below) of the Principal Field Emission Elliptic Function v. This function v is obtained by applying 
the simple- JWKB approximation [5, 6] to a Schottky-Nordheim (SN) barrier, and can be defined 
by a JWKB integral (see, for example, [7]). Although always called an "elliptic function", v is 
better described as a special kind of complete elliptic integral. An expression for v in terms of the 
complete elliptic integrals of the first and second kinds, K(m) and E(m), where m is the elliptic 
parameter as used in [8], has been known for many years, and numerical values have long been 
available [9]. 

Recently, [7, 10], an exact series expansion (see below) has been found for v, in terms of the elliptic 
variable I' defined by 

l' = [(l-m)/{l + m)f. (2) 

It follows that 

m(0 = (1 - VP)/(1 + W). (3) 

The primed symbol V was chosen, in accordance with a normal convention in elliptic-function the- 
ory, because V is a "complementary" elliptic variable in the sense that V — ► as m — > 1. Knowing 
the form of the exact expansion for v(l') has enabled us to develop [7] numerical approximation 
formulae for v(V) and dv/dl', with absolute error |e| < 7 x 10~ 10 , that substantially outperform 
earlier numerical approximations of equivalent complexity [11]. It also provides a mathematical 
explanation for the success of the simple approximation formula [7, 10] 

v(V) w l-Z' + (l/6)Z'lnZ', (4) 

which has the merits that it is exact at I' = and I' = 1, and (when assessed over the whole range 
< I' < 1) has absolute error |e| < 0.0025 and outperforms all existing formulae of equivalent 
complexity [7]. 

In CFE, the physical parameter / (the "scaled barrier field") is defined by / = F/F4, where F^ is 
the barrier field necessary to reduce to zero a barrier of zero-field height </>. The values of F and 
/ determine the strength of the tunnelling barrier. For the SN barrier, this physical parameter / 
can be identified with the mathematical parameter I', and /i in eq. (1) is given by v(V = /). 
The complete elliptic integrals of the first and second kinds may be defined [8] in terms of the 
elliptic parameter m by 

r 7 '/ 2 A(h r 71 / 2 / 

K{m)= I v = and E(m) = / J 1 - msm 2 <j)d<j). (5) 

Jo \J 1 — m sin 2 <p Jo 

Using (3), K and E can, in principle, be expressed in terms of V rather than m. The expression 
for v(V) in terms of K and E then is [7]: 

v(l') = (1 + VF) 1 / 2 [E{m{l')) - VFK(m(l'))} (6) 
where m{V) is given by equation (3). 

The exact series expansion for v(V) was originally generated [10] by the algebraic manipulation 
package MAPLE, by expanding an expression equivalent to (6). Numerical values obtained from 
MAPLE agreed with values obtained by numerical evaluation of the relevant JWKB integral, 
to better than 12 decimal places, strongly indicating that both methods were mathematically 
sound. We subsequently found algebraic formulae that reproduced the MAPLE result (see [7], 
Appendix 1). However, "manual" derivation of higher-order terms in v{V) is excessively laborious; 
also, this method does not bring out the underlying mathematics. 

Seeking a better mathematical derivation, we found [7] that v{V) is a particular solution of the 
ordinary differential equation (ODE) 

!'(l-0 5=^. (7) 
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when n = 3/16 and certain boundary conditions (see below) are satisfied. Equation (7) appears 
to be new to mathematical physics and elliptic-function theory. Its solution has not been found 
in easily available mathematical literature, but does not involve any new mathematical principles. 
Applying the method of Frobenius [12] generates a series expansion with coefficients defined by 
recurrence relations (see below). These formulae correctly generate the numerical coefficients 
already reported [10]. 

This finding has itself been reported [7] , but the full proof was too lengthy to include in earlier 
work. This paper presents the full mathematical derivation of the relevant recurrence relations. 
In part, the aim is to provide closure on the mathematics of the (simple- J WKB treatment of 
the) SN barrier; but we hope this method might also be a paradigm for finding exact solutions 
for correction factors for barriers of other shapes. It would be applicable where a method can 
be found for converting the JWKB-typc integral definition of the correction function into an 
equivalent differential equation. 

The structure of the paper is as follows. Section 2 derives boundary conditions that v must satisfy, 
Section 3 derives general solutions for (7), Section 4 uses these to derive a series expansion for v, 
and Section 5 provides discussion. 



2 Boundary conditions 



To establish the expansion for v, we need the boundary conditions that v(l') and dv/dl' satisfy at 
I' = 0. It is well known in CFE theory that v(0) = 1. However, dv/dl' becomes infinite as I' — > 0, 
so the boundary condition on dv/dl' has to take the slightly unusual form that u dv/dl' becomes 
infinite in the correct way". We develop both conditions from a result proved by Cayley in his 
1876 textbook [14]. This section is presented in the form of three lemmas. 



Lemma 1 As I' approaches zero from above, the function K(m(l')) is given by 

K(m(l')) = (3/2) In 2- (1/4) InZ' + 0(VF). 

Proof. From (3), we have 

(1 - m) = 

Formula 17.3.26 in [8] is derived from Cayley's result and states that 



(8) 



lim K(m) = In H/J(l - m) 



Hence 



lim K(m(l')) = lim In 



'l + VF 

2y/F 



lim 



(3/2) In 2- (1/4) In I' + 0(W) 



as required. We have thus proved a Cayley-type result for K as a function of I'. This result is key 
to deriving the series expansion for v. 

Lemma 2 v(0) = 1. 

Proof. Using lemma 1, we find that, in the limit of small I', the term involving K(m(l')) in 
definition (6) takes the form 



lim 



(1 + VF) 1/2 VF[(3/2) In 2 - (1/4) Inl'} 



= 0. 



When I' = 0, then m = 1, and the term in E(m(l')) in definition (6) reduces to 
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,tt/2 

E(m = 1) = / cos (j> dcj) = 1. 
Jo 



It follows that u(0) = 1. This is a well-known result in CFE theory, but for completeness we have 
given formal proof here. 



Lemma 3 



lim [dv/dl 1 - (3/16) In/'] = -(9/8) In 2 



(9) 



Proof. From [7] we have the result that 

dv 3K(m(l')) 



dl' 



1/2- 



(10) 



Using (8), we have, for small I', that 

dv ^ 3 
dF ~ ~4 



1 - 
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-ln2- -Inl' 
2 4 



The limiting form for dv/dl' in (9) follows. 



3 General solution of the ODE 

We now solve the ODE. Mathematically, it is convenient to solve the general equation, i.e. (7), 
and then put n equal to 3/16 . We use the symbol W to denote a general solution of the relevant 
equation, whether the index is taken generally as n or specifically as 3/16. 

We use the method of Frobenius [12, 13] to find a series expansion for W(V). Initially assume a 
power series solution of the form l' r X)i=o 5^'% where r and the gi are to be determined. Substi- 
tuting into (7) and equating powers of we find 



.go? 1 



(r - IT' 1 + l' r J2 + r - l)(i + r)g t - n 9i + (i + r)(i + r + l)g i+1 ] = 0. (11) 



i=0 



Since this is an identity for all powers of I', the first term has to be identically zero and there is a 
recurrence relation for the coefficient gc 



(i + r - + r) + n 

9i+i = r , w- , ~~r\ 9ii * > 0. 

(i + r)(i + r + 1) 



(12) 



Obviously, if go = then all gi (i > 1) are also zero, and W is identically equal to 0; so the go = 
case cannot contribute to the general solution. Hence, we deduce that the indicial equation is 
r(r — 1) = 0. This equation has the two roots r\ = l,r2 = 0. Since \r\ — r2 1 is an integer, two 
linearly independent solutions Wa, Wq are [12, 13]: 



W A (l') = aW 1 (l% 



(13) 



W c {l')= X [W 2 {l') + W 1 {l')\nl'], 



(14) 



where a and \ are arbitrary constants. We take a equal to unity, without loss of generality, 
but leave x to be chosen appropriately later. W\(V) and W2Q,') are the series expansions that 
correspond to the values r = l,r = 0, respectively, and are written in the forms 



(15) 



i=0 



i=0 
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where the a, and bi are coefficients to be determined. 
The general solution of (7) can thus be put in the form 



W{1') = AW A {1') + CW c (l') 



(16) 



where A and C are arbitrary constants to be determined by the boundary conditions, and the two 
independent solutions have the forms 



W A (l') = J2a i l H+1 . 



(17) 



i=0 



W c (l') = x 



^bil H + l'\nl'^a l l' i 



,i=0 



i=0 



(18) 



Note that the multiplier I'lnl' is an intrinsic part of the correct mathematical solution of (7). 
The coefficients a, have the properties of the coefficients gi when r is taken as 1 in (11). So, 
from (12): 

i(i + 1) + n 

°» 1= (i + l)(i + 2) «> (19) 

As far as the solution Wa is concerned, the value of ao can be chosen arbitrarily without loss of 
generality, because all values of ai (i > 1) are proportional to ao. 

To determine the coefficients bi, we substitute form (18) into (7), eliminate the terms in In?' by 
noting that W\ is a solution [hence, — l')d 2 Wi/d I' 2 — nW\ = 0], expand W 2 (l') in series form, 
and eventually obtain 



n(b + hi') - 2b 2 V + J2 ([{(* + + 2) + n}b l+2 - (i + 2)(i + 3)6 J+3 ] l H+2 ) + ^ cJ H = 0. (20) 

i=0 i=0 

where 

oo oo 

^c/ J = 2(Z' -l)dWi/di' = (/'-l)^(2i + l)a/\ (21) 

i=0 i=0 

Hence 

c = — a and Cj = (2i — l)aj_i — (2i + l)aj, i > 1. (22) 

In order for (20) to be an identity, we need the coefficients of l'° and l n to be zero, which requires 
that: Co + nb — 0; and c\ + nb\ — 2b 2 = 0. This yields 

b = -c /n = a /n, (23) 

b 2 = (ci + nh)/2 = (a - 3ai + n6i)/2. (24) 
More generally, the requirement that in (20) the coefficient of l n be zero yields 

{(i + l)(i + 2) + n}b i+2 - {i + 2)(i + 3)b i+ 3 + c i+ 2 = 0, i > 0, (25) 
which on re-arrangement yields 

c, + {(i - l)i + n}bj . 

= WTY) ' 1 > 2 ' (26) 

Equation (26) also gives b 2 correctly; so, using (22), we can write 

(2t - 1)^-1 - (2t + IK + {(* - l)j + ... / 97 x 

&i+i = , *>1. (27) 
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In summary, the value of 60 is fixed by (23), and the values of bi (i > 2) are determined, via (27), 
by the choice of b\. 

If Wc is to be a linearly independent solution of (7), then any multiple of Wq must also be 
a solution. We can satisfy this requirement by ensuring that each of the coefficients dj,&j is 
individually proportional to a . This is already true for all the otj and for 6 - It is sufficient that 
we further put 

61 = /too, (28) 

where /3 is a constant to be chosen; from (27), this ensures that all bi (i > 2) are proportional to 
a . 

At this point, we have derived general forms for two linearly independent solutions of (7). To 
derive specific forms for these independent solutions, we need to choose values for ao, f3 and \. In 
principle, these might be chosen in various ways. We choose them here in such a way as to make 
the boundary conditions at /' = straightforward to apply. 

From equations (16) to (18), and (23), the form for W{1'), up to terms involving I', is: 

W(l') = Caox/n + {Aa + C X pa )l' + C X a l' In J' + ... , (29) 
and the corresponding form for dW /dl' is 

dW/dl' = Aa + C/3xa + C X a + C X a In I' + . . . (30) 

From (29), when V = 0: 

C = (n/ X )(W /a ). (31) 

where Wq=W(0). Clearly, a convenient way to deal with (30) is to choose (3=—l. Equations (29) 
and (30) are now reduced to 

W(l') = W + [Aa - nW ]l' + nW l' In I' + . . . (32) 

dW/dl' = Aa + nW In l' + . . . (33) 

Since these equations contain ao only in the combination Aao, we can take ao — 1 without loss of 
generality; then b = l/n and b\ = —1. Values of a, (i > 1) and bi (i > 2) can then be obtained 
from the recurrence formulae above. 
At this point, we have the situation that Wa(0) = and 

Wc(0) = X bo = x/n (34) 

It is convenient to have Wc(0) = 1, so we choose x — n > which also simplifies (31) to 

C = W(0)/a o = W(0). (35) 

This completes the definition of two convenient, linearly independent solutions of ODE (7), and 
enables the general solution to be written in the form 

00 

W{1') = l + Y^{Aa i + Cnb i+1 +Cna i \nl')l' i+1 , (36) 

i=0 

where a = l,b\ = —1, and the other coefficients are given by the recurrence relations (19) 
and (25). Note that this form (and the value derived below for C) differ slightly from a form 
proposed earlier [15], because we have defined Wc in a slightly different (more convenient) way 
here. 
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4 Series expansion for function v 

At this point, we revert to discussion of the ODE of index 3/16, and derive the expansion for v(l'). 
For n = 3/16, the first few terms of the two independent solutions are 

^ (0 = r + i,« + T | ji « + ..., (3 7, 

= 1 - r/ + iii'' 2 + g&r + -) ''' n ''- <38 > 

Using the boundary condition Wo = v(0) = 1, (35) yields C = 1. For the second boundary 
condition we use (as noted earlier) the slightly unusual requirement that, as /' tends to zero, 
the derivative dW/dl' "tends to infinity in the correct way". Comparing (9) and (33) yields the 
following identity, valid in the limit of small /'-values: 

(3/16) In/' - (9/8) In 2 = nW In/' + Aa (39) 
Since n = 3/16 and Wo = 1, the terms in In/' cancel. So, since do = 1: 

A = -(9/8)ln2. (40) 

So the first few terms of the expansion for v(V) are: 

„n - /9, n 3 \ / 27 , „ 51 \ „, / 315 , „ 177 \ „, 

w =1- - In 2 H /'- In 2 In 2 Z' 3 + . . . 

y ' \8 16 J V256 1024/ V8192 8192/ 

V16 512 16384 J K ' 

This series has the form reported earlier [7, 10], found by using MAPLE. Evaluating the coefficients 
to 5 decimal places yields 

v(l') = l-0.96729/'-0.02330/' 2 -0.00505/' 3 +. . .+/'ln/'(0.18750+0.01758/'+0.00641/' 2 +. . .). (42) 

By dividing the first expression by 1 — /', this can be put into a form that is exact at both /' = 
and /' = 1, even when the two series expansions below are truncated: 

v(l') = (l-/')(l + 0.03271/' + 0.00941/' 2 + . ..) + /' ln/'(0.18750 + 0.01758/' + 0.00641/' 2 + . . .). (43) 

Both of the series in (43) have good convergence properties in the range < /' < 1, so this form 
is useful for numerical approximations [7]. 

5 Discussion 

For most of the last 50 years, the correction function v has been expressed in CFE theory as a 
function of the Nordheim parameter y — \fV [4], rather than /'. In [10] it was argued that, because 
the exact series expansion (41) contains no terms in VI 7 , it is mathematically more natural to use 
/' as the independent variable. The derivation here, using the method of Frobenius, confirms this 
formally. Terms in In/' are an intrinsic part of the correct mathematical solution, but fractional 
powers of /' are not necessary. 

Approximate formulae such as (4) cannot be derived by simple Taylor expansion methods. It is 
implicit in the mathematical analysis here that simple Taylor expansion methods cannot generate 
good approximate formulae valid for the whole range < /' < 1. This conclusion may also be 
applicable to tunnelling barriers of other shapes. 

Equation (7), the defining ODE, can be transformed to be in terms of the variable y. This yields 

2 ,d 2 W l-y 2 dW 3 Tir n 

(l-2/)-ro- + — ; t^ = 0. 44) 

dy 2 y dy A v ' 
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Obviously, this equation is more complex than (7). This is a further good mathematical reason 
for choosing to use the variable I', rather than y, in CFE theory. As discussed elsewhere [7], an 
important implication of the role of I' in the mathematics is that the natural variable to use in 
physical discussions of CFE theory is the scaled barrier field /. 

The algebraic manipulation package MAPLE has played a crucial role in stimulating this work, 
because our use of MAPLE [10] drew attention to the existence and form of the series expansion 
for v(l') and provided a result to aim for. Derivation of this series by inserting series expansions for 
K{1') and E{V) into (6) proved excessively laborious if performed by hand [7], even when we had 
found the Cayley forms [14] for the series expansions of K(m) and E(m), (which, interestingly, 
are not reported in several of the elliptic- function handbooks consulted). So we looked for the 
ODE that v(l') obeyed and found (7) [7]. 

Luckily, the ODE (7) has a very simple form. Nevertheless, the derivation of solutions is quite 
a lengthy mathematical process. Wylie [12] (p. 322) describes the solution of equations of this 
general kind as "straightforward but very tedious" ; we confirm his description. Much of Section 3 
here could, in principle, have been written in the late 1800s, or at any subsequent time, but there 
has been little incentive to do so. 

Overall, it is not entirely surprising that getting a series expansion/definition for v in place has 
taken nearly 80 years, measured from the original (incorrect) attempt by Nordheim [4] to derive an 
exponent correction function for the SN barrier. Essentially, it is the relatively recent introduction 
of reliable computer algebra packages that has made the discovery of this expansion (and its 
mathematical origin) feasible, rather than so difficult as to be unlikely to happen. The problem 
for CFE and the SN barrier (even when using the simple- JWKB approximation) has been that the 
necessary mathematical functions have not been solidly in place when physicists needed them, in 
the 1920s and 1950s; one might contrast this with the wave-mechanical analysis of the hydrogen 
atom in the 1910s, which needed mathematical functions that had been well defined for very many 
years. 

As indicated in the introduction, our hope is that the derivation in this paper may be able to serve 
as a paradigm for the treatment of other barrier models, particularly for realistic models for the 
potential-energy variation above sharply curved emitters. The keys, in each case, will be to find 
an ODE that the tunnclling-exponent correction function satisfies, and a suitable formulation of 
the boundary conditions. 
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